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A quantum algorithm can find a target item in a database faster than a classical algorithm. One 
can trade accuracy for speed and find a part of the database (a block) containing the target item 
even faster, this is partial search. We consider different partial search algorithms and suggest the 
optimal one. Efficiency of an algorithm is measured by number of queries to the oracle. 



I. INTRODUCTION 

Database search has many applications and is used widely. Grover discovered a quantum algorithm that searches 
faster than a classical algorithm [1]. It consists of repetition of the Grover iteration G\, which operates on the 
computational quantum states. The number of repetitions is: 

ifull = jS* (!) 

for a database with large number of entries N . After jf u \\ the algorithm finds the target item. For more details, see 
also [2-4]. Below we shall call G\ a global iteration. 

Sometimes it is sufficient to find an approximate location of the target item. A partial search considers the following 
problem: a database is separated into K blocks, of a size b = N/K. We want to find a block with the target item, not 
the target item itself. Such partial search was first introduced by Mark Hciligman in [5] , as a part of algorithm for list 
matching. We can think of partial search in following terms: an exact address of the target item is given by a sequence 
of n bites , but we want to find only first k bites (k < n). Fast quantum algorithm for a partial search was found by 
Grover and Radhakrishnan in [6]. They showed that classical partial search takes ~ (N — b) queries, but quantum 
algorithm takes only ~ (V~N — coefiS/fe) queries. It uses several global iterations Gj 1 and then several local iteration 
&2 , see (8). Local searches are Grover iterations [searches] in each individual block made in each block separately in 
parallel. Grover- Radhakrishnan algorithm was improved and simplified in [7]. The number of queries to the oracle 
in this algorithm was minimized by in [8], the coeff was maximized. Below we shall explain minimized version of 
Grover-Radhakrishnan algorithm. We shall call it GRK algorithm. In this paper we consider three other versions of 
partial search algorithm. They use different sequences of global and local searches: local-global, global-local-global 
and local-global-local. We prove that GRK version still uses minimal number of queries to the oracle. We conjecture 
that GRK algorithm is optimal among all partial search algorithms, which consist of arbitrary sequence of local and 
global searches. 

The plan of the paper is as follows. In the next section we remind the Grover algorithm. After this we formulate 
minimized version of Grover-Radhakrishnan algorithm [GRK algorithm] . In the rest of the paper we consider other 
partial search algorithms. We arrive at the conclusion that GRK uses minimal number of queries comparing to other 
algorithms. 

II. PARTIAL SEARCH 
A. Global Iterations 

First let us remind the full Grover search. We shall consider a database with one target item. The aim of the 
Grover algorithm is to identify a target state |i) among an unordered set of N states. This is achieved by repeating 
global iteration which is defined in terms of two operators. The first changes the sign of the target state \t) only: 

I t =I-2\t)(t\, (t\t) = l, (2) 

where / is the identity operator. The second operator, 

4 =/-2| Sl )( Sl |, (3) 
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changes the sign of the uniform superposition of all basis states \s\), 



N-l 



x=0 



The global iteration is defined as a unitary operator 

G\ = —I Sl It- 

We shall use eigenvectors of G\: 

G 1 |^ 1 ± ) = A±|Vr t ), Xt = exp[±2i0 1 ], \tf) = l=\t)±±= 
They were found in [4], where the angle 6\ is defined by 



JV-l 



Wt 



sin 2 6»i 



N' 



(4) 



(5) 



(6) 



(7) 



B. GRK Algorithm for Partial Search 

The first version of partial search was found in [6]. The algorithm uses j\ global iteration and ji local iterations. 
LOCAL ITERATIONS are Grover iterations for each block: 



G2 — —I S2 It- 

It is given by (2), but I S2 is different. In one block it acts as: 

1 ^\block = i \uock- 2 \ S ^^ l*2> 



(8) 



(9) 



one block 



In the whole database I S2 is the direct sum of (9) with respect to all blocks. Both relevant eigenvectors of Gi were 
found by in [4]: 



G 2 |V 2 ± )=A±|V 2 ± ), Xf = exp[±2i6 2 ], \4,f) = —\t)±—\ntt) 
Here the |ntt) is a normalized sum of all non-target items in the target block: 

|ntt) = -=L= Yl (ntt|ntt) = l. 



target block 



We shall need an angle 62 given by 



■ 2n K 1 



The partial search algorithm of [6] creates a vector 

\d) = G 1 Gi 1 G{°\s 1 ). 



(10) 



(11) 



(12) 



(13) 



see 1 . In the state \d) the amplitudes of all items in non-target blocks are zero. Notice that this algorithm uses 
global-local sequence of searches. We consider large blocks b = N/K — > 00. The number of blocks K is an important 
parameter. We shall replace it with 



1 We use a modification of [6], suggested in [8] 
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1 7T , 

sin7=— 7=, < 7 < - (14) 



The optimal version of this algorithm was find in [8]. It can be described by the following equations: 

n \ sin 7 cos 27 . n , cos 27 , „ N 

cos 2 ji 2 ) = ' . tan(2j 6 l i) = f (15) 

V 21 cos7sin2 7 ' v J0 11 (5^7)^3-4(8^7)2 V 7 

Partial search is faster then full search (1) by ~ y/b, the coefficient in front of \fb is explicitly calculated in [8] and 
studied as a function of number of blocks. 



C. Notation and Setup for General Partial Search Algorithm 

Let us introduce a unite vector: 

IU)= Jb(K-l) ^ H«> = 1- (16) 

all items in all 
non-target blocks 

We shall use a three dimensional space. The orthonormal basis is formed by the target item |t), sum of all non-target 
items in the target block \ntt), defined in (11) and \u). All the state vectors involved in present quantum search 
problem can be written in this basis as 

\V>=(a,b,c) T (17) 

meaning 

\V >= a\t > +b\ntt > +c\u > (18) 
For example, the global uniform state which is used as initial state of searching is given by 

|si) = (sin 7 sin 62, sin 7 cos 62, cos 7) T (19) 

and the local uniform state is 

1*2) = (sin 6» 2 , cos 6» 2 , 0) T (20) 

The algorithms which we consider in this paper can be represented as matrices in this linear space. For example j'2 
repetitions of the local iteration (8) is: 

/ cos(2j 2 # 2 ) sin(2j 2 # 2 ) 0\ 
G J 2 2 = -sin(2j 2 2 ) cos(2j 2 2 ) (21) 
V ij 

The ordering of eigenvectors is \t), \ntt) and \u). The matrix has three eigenvectors: 

G¥ \vf) = cM±^9 2j2 )\vf), \v° 2 ) - \v° 2 ) (22) 

The eigenvectors can be represented as: 

^=7f(t ! )' !*-(;)- (23) 

Now let us turn our attention to global iterations (5), j\ repetitions of the global iterations can be represented as 

cos(2ji#i), sin(2ji^i) sin7, sin(2ji#i) cos 7 

Gf = [ -sin(2ji6»i)sin7, cos 2 7 + cos(2j 1 1 ) sin 2 7, sin7COS7 + cos(2ji6»i)) ] (24) 

- sin(2ji#i) COS7, sin7COS7 ((— l) 3l+1 + cos(2ji#i)) , (— l)- 71 sin 2 7 + cos(2ji#i) cos 2 7 
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This is a simplified asymptotic expression valid in the limit of large blocks b — > oo. We used (14). The matrix has 
three eigenvectors: 

Gj 1 ^) = ex p{±2i0 1 j 1 )\vt), Gj>?) = (25) 
The eigenvectors can be represented as: 



\vt) = 4s ( ±* sin 7 ] . \v°) = ( COS7 ] . (26) 
±i cos 7 / \ — sin 7 



V2 



D. General Problem of Partial Search 



The ultimate goal of partial search is to start with the uniform state \s\ > and locate the target block, and obviously 
GRK is not the only means to achieve the goal. Based on the two types of queries, global and local iterations, we 
naturally generalize GRK into a wide set of partial search algorithms by alternate use of the two iterations: 

G(j k ,j k -u ■ ■ ■, j 2 , ji, jo) = G^G^Gi^G^' 3 ■ ■ ■ G 1 i 'G 2 h G 1 io (27) 
which fulfills the following condition 

< u\G{ju,3k-\, - ■ -,j2,ji,jo) \si >= (28) 

This means that amplitude of each item in non-target block is zero. Here all ji are non-negative integers. The total 
number of queries for these algorithms is given by S — Yli=ofri an d we should try to minimize S to find the optimal 
one. 

Now let's consider various sequences. A few discussions can be made here: 

1) The last queries in the sequence should be the global iterations. Note that local iteration (21) doesn't do anything 
on \u > but only rotates state component inside the target block, so if an algorithm makes computational state in the 
target block after the last local queries, it can simply waive the last local queries since the state must already rest in the 

target block before those unnecessary last local iterations. That is, if < u\ G^^ G\ h G^ 1 ■ • -G\ 2 G^ 1 G\° \ s \ > = 

then we must also have < u\ G\ h G^* 1 ■ ■ ■ G\ 2 gV* Gi° |si >= 0. 

2) The first queries could be either local or global iterations. Sequences starting with the global (jo > 0) include: 

Gi (Global), Gi G 2 Gi (Global-Local-Global), and so on. The simplest one, with only global queries involved, 
gives nothing but the original Grover's full search (1), which saves no steps but precisely locates the target item. The 
next simplest one, Global-Local-Global, will be studied in later section. Note that the up-to-now established optimal 
partial search GRK algorithm also falls into this category, with j 2 = 1 and jo, ji specified by (15). 

3) For those starting with local queries (jo = and ji > 0), we have: d (Local-Global), Gi Gfd G^ 1 
(Local-Global-Local-Global), and so on. We will later fully consider the Local-Global algorithm and also discuss one 
particular case of Local-Global-Local-Global which has j\ to be one, namely with only one global query applied after 
the last local iterations (we call this as Local-Global-Local sequence). 

4) The argument based on the optimality of Grover's full search puts lower limit of the partial search steps S > 
l\fN - \\fb. It was shown in [6] that S — jVn — R^/b. Here R is a constant independent of b. The optimal 
scheme we are seeking should maximize R. Also according to this thought, we expect the total number of global 
queries should scale asymptotically as jVn — r/Vb, since in partial search we should be faster than full search j 
but save steps only of order \fb. As for total number of local queries it should scale like aVb, since irVb times local 

iteration will have rotated the state vector a whole lap on the target block \t > \ntt > plane. It follows from 

these consideration that R~rj — a. Though the numbers of local and global iterations ji are integer numbers and in 
principle not continuous variables, in the large block limit we can reasonably treat them as quasi-continuous and use 
them as function arguments. 

5) R is a function of block number K or instead the parameter 7 = arcsin(-^=) with < 7 < 7r/4. Another limit 
we will consider is the large K limit, K >> 1 or cquivalcntly 7^0. In this limit the mathematical formulae will be 
significantly simplified, which will be very helpful to analytical efforts as will be seen in later sections. The cases with 
not so large K could be complemented by direct numerical verification. 
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III. LOCAL-GLOBAL SEQUENCE OF SEARCHES 

The Local-Global sequence of searches is the simplest partial search scheme. For such a sequence G(j 2 ,ji) with 
first ji local iterations applied and then j 2 global, we have the final computational state to be: 

G 1 h G 2 h \s 1 >={A,B,C) T (29) 

Here A,B,C are the coefficients of components \t >, \ntt >, \u > respectively (see (18)), which are given by 

A = sin-f coa(2j 2 9i) sin(2ji8 2 ) + sin 2 7sin(2j 2 #i) cos(2j 1 2 ) + cos 2 7 sin(2j 2 ^i) 
B = - sin 2 7sin(2j 2 ^i) sin(2j 1 6» 2 ) + (-) h cos 2 7sin7Cos(2j 1 6» 2 ) 

+ sin 3 7cos(2j 2 6» 1 )cos(2j 1 2 ) + shi7 cos 2 7 (cos(2j 2 6> 1 - (-)- 72 )) 
C = — sin 7 cos 7 sin(2j 2 #i) sin(2ji# 2 ) + sin 2 7 cos 7 (cos(2j 2 #i) — (— ) J2 ) cos(2ji# 2 ) 

+ (-) h sin 2 7 cos 7 + cos 3 cos(2j 2 6>i) (30) 

Then to accomplish the partial search we should have the constraint equation 

< u\G 1 n G 2 1 \s 1 >=C = (31) 

which means that the amplitude of every item in non-target block vanish (28). Since there are only one set of local 
and one set of global iterations, we apply the scaling as discussed before and introduce 

h = aVb , j 2 = ^VN - r)Vb (32) 



We have the limitation < a < n and rj < jVK. Now we can rewritten the constraint equation as following: 

( 277 2?7 \ 

cos7- I [sin 2 7Cos(2a) + cos 2 7] sin(^^) — sin7sin(2a) cos(^^) + (— )- 72 sin 2 7[1 — cos(2a)] I = (33) 

Note that K > 2 so cos 7 > 0, hence it drops out in the above equation and we have 

F(a,r]) = [sin 2 7 cos(2a) + cos 2 7] sin(-p=) - sin7sin(2a) cos(^4) + {-) 32 sin 2 7[l - cos(2a)] =0 (34) 



Remember we are interested in the fastest algorithm (using the fewest queries), especially those using less steps 
than Grover's full search. Now we have the total number of queries to be S = jVN — (rj — a)Vb, thus to minimize S 
we should maximize R = rj — a under the constraint (34). The point here is that from the constraint we can consider 
77 as a function of a, thus R is also a function of a, to which it should be optimized. We then have ^ = — ^ + 1 = 0, 
which leads to 

2n 

sin(27)[cos 2 7 (1 - cos(2a) cos(-p=)) + {-) n sin 7 sin(2a)] = (35) 

v K 

Again since < 2 7 < n/2, we have 

sin(2 7 ) > (36) 

so we simplify the above equation as 

2ri 

H(a,rj) = 2 sin a [cos 2 7 sin a cos (— + (— Y 2 sin 7 cosa] = (37) 

yK 



There is an apparent solution for (34) and (37), namely sin a = sin(-^=L) = 0. This gives a — and rj — 1-tv-\/K, I = 
0, — 1, —2, • • •, but to maximize R = rj — a we should have a = r\ = 0, which again recovers the Grover's full search 
solution and is a trivial one for partial search. In the following we only consider nontrivial solutions with sin a 7^ 0. 

Let's first look at the large K or small 7 limit. In leading order, the two equations (34) and (37) are reduced to be 

F(a, r,) = sin(-^L) - 7 sin(2a) cos(-^=) + 7 2 (-l) j2 2 sin 2 a = (38) 
\l K v K 
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H(a, V ) 



2n 

sina[sinacos(^^) + (— iy 2 -fcosa] = 
v K 



(39) 



To satisfy (38), the first term sin(-^|) must be at least as small as <~ 7, so in leading order we must have cos(-^=) <~ 1. 

This, combined with (39), requires also sin a must be at least as small as ~ 7 (we don't consider sin a — as mentioned 
before), which means cos a ~ 1. So finally we reduce (39) in leading order to be 



sin a + (— l) j2 7 = 



(40) 



Remember we have < a < ir thus sin a > 0, so the above equation has solution ONLY for ODD values of ji- For 
odd ]2 the solution is sin a = 7, this yields two solutions a = 7 and a = ir — 7. By substituting sin a = 7 back into 
(38) we get in leading order sin -^L = 2 7 2 , which means -^L = 2 7 2 + 2Ztt or -^L = -2 7 2 + (21 - 1)tt, I = 0,-1, -2, • • •. 
But again remember our goal is to maximize R = r\ — a, so we adopt the optimal Local-Global solution a — 7 and 
i] = 7 in large K limit. This, however, seems giving no speedup compared with full search since R = r\ — a ~ 0. To 
clarify this, we need go to higher order of equations (38) and (39). By properly including corrections up to <~ 7 s we 
can find the solution to be 



: 7 + - 7 3 + o( 7 4 ) 



V = 7 + g7 3 + o(7 4 ) 



Now we see in large K limit, the Local-Global sequences of search do achieve speedup 

R=\l* + otf) 

which is very small but still nonzero. 

Let's now turn back to finite values of K. From (37) we have 



cos(^=) = 

Wk' 



(-) 



j'2+1 



sin 7 cos a 



cos z 7 sin a 



Combine (43) and (34) together we get 



, 2n . sin 7 cos a 

cos(^=) = 5 : 

sJK cos^ 7 sin a 



. 2r) 2 sin 2 7(cos 2 7 sin 2 a + cos 2 a) 

Slllf — • — ) = 

\[~K cos 2 7[sin 2 7 cos(2a) + cos 2 7] 



(41) 



(42) 



(43) 



(44) 



By requiring cos 2 (-^=) + sin 2 (-^=) = 1 we obtain the equation determining the value of a and from a. we can obtain r\. 
Generally it is hard to analyze the problem analytically, so we proceed to numerical solutions , see the figures Fig. 1,2. 
In these figures we plot R as a function of a with different K for both even and odd values of J2, at each point rj is 
determined from the constraint (34) . For the odd ji case there is always a positive maximum which is faster than full 
search, and as K approaches very large values, the maximum of the e v.s. a curve moves gradually toward the origin 
which stands for the full search solution. But for even j 2 , R is always negative, hence slower than full search. 




-0.2 



0.5 



1 . 5 
a 



2.5 




FIG. 1. Dependence of R on a with odd j'2 in Local-Global sequence, the curves from top to bottom are for 
K = 5,6,7,8,9,18,36,72,144,200 respectively. For each a the value of ij is solved from (34). Right panel is the amplifi- 
cation of the area near origin in Left panel. 
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0.5 1 1.5 2 2.5 3 

a 

FIG. 2. Dependence of R on a with even j2 in Local-Global sequence, the curves from top to bottom are for 
K = 5, 6, 7, 8, 9, 18, 36, 72, 144, 200 respectively. For each a the value of 77 is solved from (34). 

An important comparison is to be made between this Local-Global algorithm and the established GRK one. Our 
numeric results show that the odd j 2 Local-Global searches can get R w 0.3 with K = 2, 3 and get < R < 0.09 
for K > 4. For GRK, however, it can achieve R > 0.32 for all K > 2. So the GRK is much faster than the present 
Local-Global searches. 

To sum up with the Local-Global sequence of searches, we find it can be faster than full search with R\/b speedup 
for block number K up to several tens, but it is always slower than the GRK optimized partial search. 



IV. GLOBAL-LOCAL-GLOBAL SEQUENCE OF SEARCHES 

In this section we shall consider another algorithm for partial search, the Global-Local-Global sequence of searches, 
which starts with (4), first applies jo global iterations (5), then j\ local iterations (8) and finally j 2 global iterations: 

G^G^Gi 30 ^). (45) 
Again amplitudes of all items in non-target blocks should vanish at the end of algorithm: 

(u\G 1 j2 G 2 il G 1 i °\s 1 )=0. (46) 
Let us present this equation in an explicit form. After first set of global iterations the state of the database will be: 

s - / cos(2j o 0i) \ / sin(2j o 0i) 
Gi J °|«i) = —J- -sin(2j tfi)sin 7 + S in 7 cos(2jo0i) | . (47) 
y — sin(2j #i) COS7 J y cos7cos(2jo#i) 

After local iterations the state of the database is: 

sin(2j o 0i) cos(2ji# 2 ) + sin7 cos(2j #i) sin(2ji6> 2 ) 
!ji02) + sin 7 cos 
cos7cos(2jo#i) 



G 2 n Gi°\s 1 ) = I -sin(2joei)sin(2jie 2 )+sin7cos(2jofli)cos(2jie 2 ) I. (48) 



Here we neglected ~ 1/y/b terms (namely taking large block limit). After next set of global searches we have to 
calculate only third component (coefficient of \u)) of the vector: 

0=(u\G 1 j2 G 2 jl G 1 jo \s 1 ) (49) 
= cos(2ji(9 2 ) {-cos7sin(2j 2 <9i)sin(2.)o6»i) + sin 2 7cos(2j 6'i) cos 7 [ (-1) J2+1 + cos(2j 2 6>i) ]} + 
sin(2ji^2) {— sin7COS 7 sin(2j 2 6'i) cos(2j #i) — sin7COS7sin(2j 6'i)[ (— 1) J2+1 + cos(2j 2 ^i) ]} + 
+ cos7cos(2j #i)[(— 1) J2 sin 2 7 + cos 2 7 cos(2j 2 6q)] 
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This constraint equation guarantees the amplitudes of all items in all non-target blocks vanish to successfully complete 
the partial search. 

First let us check the case j\ = 0, no local searches. In this case the constraint equation can be reduced to 
cos[26>i(jo +32)] = 0. This is just the full search: we use jo + 32 = f \/N global iterations to find the target item. 

Now let us consider more general and complicated case. We expect the following scaling, namely the global iterations 
jo = § VN — TfVb, 32 = fiVb and total local iterations ji = a\fb with r\ < j\/K, < (5 < \\f~K and < a < it. It 
should be remembered that our purpose is to minimize jo + ji + 32 = § V^V — RVb, R = r\ — (3 — a, or equivalently 
maximize R. Our strategy is similar to that used for analyzing Local-Global sequence: study the large K limit 
analytically while deal with finite K case numerically. 

In the large K or small 7 limit, we can take the leading order of (49) and simplify our constraint to a much simpler 
form : 

r] = /3 cos(2a) + 1 - sin(2a) (50) 

From this equation, there are two possibilities: 

1) 32 is even, thus r\ — /3cos(2a). In this case the problem simplifies into maximizing R = /3(cos(2a) — 1) — a with 
< < and < a < ir. Note that cos(2a) — 1 < so the maximum of R must occur at rj = (3 = a = 0, which 
is again the trivial full search solution (1). 

2) 32 is odd, thus rj = /3cos(2a) + sin(2a). We then have to maximize R = f3(cos(2a) — 1) + sin(2a) — a with 
< P < jVK and < a < it. The solution is 

7?= ^ , p = , a- tt/6, (51) 

which has achieved R^/b ss 0.3424V5 speedup with respect to full search. This nontrivial optimal solution is exactly 
the GRK algorithm, with vanishingly small odd ]2 namely 32 = 1- 

Now let us discuss finite values of K. We calculated dependence of R on /3 by determining a from constraint 
equations and optimizing 77 numerically at each value of (3. We done this numerically for up to 200 blocks. The 
dependence of R on (5 is monotonous. Corresponding figures took too much memory, so we withdraw them from the 
paper. For even j\ searches we consider are slower than full Grover search, while odd 32 these searches are faster than 
full search. The optimal searches always occur with odd and vanishing j 2 (the GRK case), being about 0.34-^ faster 
than full search. The large K case also confirm our analysis above. Our results here also confirm that the optimum 

partial search is GRK [8], it has the form GiG^ 1 G\° ■ 

To conclude, we have found that Global-Local-Global sequence of searches can be much faster than full search (1) 
with odd times global iterations applied in the end, and the optimal search is determined to be the GRK algorithm 
for all values of K. 



V. LOCAL-GLOBAL-LOCAL SEQUENCE OF SEARCHES 

In this section we shall discuss another algorithm of partial search. It belongs to the category of Local-Global-Local- 
Global sequence but with only one global query applied in the end, which we call as Local-Global-Local sequence. We 
start with (4) first apply ji local iterations, then j 2 global iterations (8) and then j 3 local iterations, and eventually 
a single global query: 

GiG 2 j3 Gi 2 G 2 n |si). (52) 

One particular advantage of this type of sequence is that when ji tends to zero our Local-Global-Local will degenerate 
to the GRK type algorithm. Amplitudes of all items in non-target blocks should vanish at the end: 

(u|GiG 2 J3 Gi J2 G 2 Jl |si) =0. (53) 

As before, we introduce the scaling of iteration numbers to be ji = aVb, 73 = SVb, and j 2 = j\^N — r/Vb. The total 
number of queries will be S = f VN - RVb with R = rj — a — 8 which we want minimize under constraint (53). The 
explicit form of the above constraint equation is as follows: 
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= X ■ sin(2a) + Y ■ [cos(2a) - l] + Z 

2ri 2n 2n 

X = sin7sin(27) sin(2(5) sin(— =) + sin 2 7sin(27) cos(2<5) cos( — =) + sin 7 cos 7 cos (27) cos(^^) 

v K v K v K 

2r\ 



(54) 



2?7 



2r] 



Y = sin 7sin(27) sin(25) cos(— — sin 7sin(27) cos(25) sin(— r=) — sin 7COS7 

VK v K 

2ti 2ti 2tj 

Z = sin(27) sin(2<5) cos(^^) — sin7sin(27) cos(2i5) sin(^^) — cos7Cos(27) sin(^^) 
v K v K v K 



cos(27)sin(- 7 ^) + (-p 



Though the above equation is complicated, it is very easy to solve numerically. Also taking large K limit can 
significantly simplify it. So as we did before, we analytically study the large K limit, complemented by numerical 
results from very small K to very large K. 

By taking large K or small 7 limit, we obtain from leading order of (54) the following: 



sin(2a) + 2 sin(25) - 2rj = 
We then have r] — \ sin(2a) + sin(25) and hence the total number of queries to be 



(55) 



R = - sin(2o:) - a + sin(25) - S 



(56) 



By requiring ^ = = we get the solution maximizing R 



a = 



6 ' 



n : 



V3 

2 



(57) 



This, with zero j\, again recovers the GRK optimized solution. So in large K limit we see Local-Global-Local sequence 
is no faster than GRK algorithm. Different from Local-Global and Global-Local-Global, here in Local-Global-Local 
we notice that the odd j'2 and even j 2 converge to each other when approaching the optimal solution and the oscillation 
terms in (54) with factor (— Y 2 disappear. 



-0 . 5 



-2 . 5 




. 3 




FIG. 3. Dependence of R on a with odd ji in Local-Global-Local sequence, the curves from top to bottom are for 
K = 5, 6, 7, 8, 9, 18, 36, 72, 144, 200 respectively. For each a the value of r/ is solved from (54) with 6 optimized numerically. 
Right panel is the amplification of the area near origin in Left panel. 
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FIG. 4. Dependence of R on a with even j2 in Local-Global-Local sequence, the curves from top to bottom are for 
K — 5,6, 7, 8, 9, 18, 36, 72, 144, 200 respectively. For each a the value of r\ is solved from (54) with 8 optimized numerically. 

For general values of K, we conduct the numerical method and show the results in figures Fig. 34, which plot R as 
a function of a with -q solved from (54) and 5 optimized numerically. As can be seen, the optimal solutions always 
occur with a — which goes back to GRK case. Also we note that even and odd j% give same results around optimal 
point. 

So in this section we have established that the Local-Global-Local sequence of searches can be much faster than 
full search, but is no faster than GRK algorithm. In the appendix an alternative approach for Local-Global-Local 
sequence based on a conjecture about cancellation of oscillation terms in (54) is briefly described, which though is not 
directly relevant here but arrives at similar result and may shed light for future exploration of even more complicated 
sequences. 

VI. SUMMARY 

We considered different partial search algorithms, which consists of a sequence of local and global searches. We 
introduced a general framework for studying partial quantum search algorithms and classified various possible se- 
quences. Particularly, we studied the Local-Global, Global-Local-Global as well as Local-Global-Local sequences of 
searches by combining numerical study for wide range values of K and analytical results for large K limit. All these 
algorithms achieve \fb speedup compared to the Grover's full quantum search. GRK algorithm [8] is the fastest among 
partial search algorithms, which we considered. 
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APPENDIX: APPENDIX: REMARKS ON LOCAL-GLOBAL-LOCAL SEQUENCE 

Let us make two remarks: 

1) Let's start from the explicit form of constraint equation (54). We minimizing j\ + ]2 + js- It is interesting to 
note that at the minimum the oscillation terms cancel. The coefficient at (— l)- 72 vanishes because of GRK equation: 

cos7sin27cos(202j3) =sin7Cos27 (Al) 

This is exactly the first equation of (15). We also can write it in the form: 

/•«•/>•< 1 - tan2 7 ■ m ■ n i V^3 - 4 sin 2 7 
cos(2j 3 2 ) = t L , sm(2j 3 # 2 ) = ~ 2 1 (A2) 

A A COS '"V 
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2) Minimum number of iterations for Local-Global-Local sequences corresponds to ji = case and the algorithm 
is reduced back to GRK version of partial search, see small ji increase ji + ji very little since d{j2 + ji)/dji = and 
d 2 (.72+Ji)Mi 2 = 0at 3l =0. 
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